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Abstract 

A connection between integrability properties and general statistical prop- 
erties of the spectra of symmetric transfer matrices of the asymmetric eight- 
vertex model is studied using random matrix theory (eigenvalue spacing dis- 
tribution and spectral rigidity). For Yang-Baxter integrable cases, including 
free-fermion solutions, we have found a Poissonian behavior, whereas level 
repulsion close to the Wigner distribution is found for non-integrable models. 
For the asymmetric eight-vertex model, however, the level repulsion can also 
disappear and the Poisson distribution be recovered on (non Yang-Baxter in- 
tegrable) algebraic varieties, the so-called disorder varieties. We also present 
an infinite set of algebraic varieties which are stable under the action of an 
infinite discrete symmetry group of the parameter space. These varieties are 
possible loci for free parafermions. Using our numerical criterion we have 
tested the generic calculability of the model on these algebraic varieties. 

PACS numbers: 05.50.+q, 05.20.-y, 05.45+b 
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I. INTRODUCTION 



Since the early work of Wigner |]J random matrix theory (RMT) has been applied with 
success in many domains of physics 0. Initially developed to serve for nuclear physics, 
RMT proves itself to provide an adequate description to any situation implying chaos. It 
has been found that the spectra of many quantum systems is very close to one of four 
archetypal situations described by four statistical ensembles. For the few integrable models 
this is the ensemble of diagonal random matrices, while for non-integrable systems this can 
be the Gaussian Orthogonal Ensemble (GOE), the Gaussian Unitary Ensemble (GUE), or 
the Gaussian Symplectic Ensemble (GSE), depending on the symmetries of the model under 
consideration. 

In the last years several quantum spin Hamiltonians have been investigated from this 
point of view. It has been found that ID systems for which the Bethe ansatz applies have 
a level spacing distribution close to a Poissonian (exponential) distribution, P(s) = exp(— s), 
whereas if the Bethe ansatz does not apply, the level spacing distribution is described by 
the Wigner surmise for the Gaussian orthogonal ensemble (GOE): 

P( S ) = ^exp(-vr S 2 /4). (1) 

Similar results have been found for 2D quantum spin systems Other statistical prop- 

erties have also been analyzed, showing that the description of the spectrum of the quantum 
spin system by a statistical ensemble is valid not only for the level spacings but also for 
quantities involving more than two eigenvalues. 

In a recent letter || we proposed the extension of random matrix theory analysis to 
models of classical statistical mechanics (vertex and spin models), studying the transfer 
matrix of the eight- vertex model as an example. The underlying idea is that, if there actually 
exists a close relation between integrability and the Poissonian character of the distribution, 
it could be better understood in a framework which makes Yang-Baxter integrability and 
its key structures (commutation of transfer matrices depending on spectral parameters) 
crystal clear: one wants to switch from quantum Hamiltonian framework to transfer matrix 
framework. We now present the complete results of our study of transfer matrices and a 
detailed description of the numerical method. This work is split into two papers: the first 
one describes the numerical methods and the results on the eight-vertex model, the second 
one treats the case of discrete spin models with the example of the Ising model in two and 
three dimensions and the standard Potts model with three states. 

We will analyze a possible connection between statistical properties of the entire spectrum 
of the model's transfer matrix and the Yang-Baxter integrability. A priori, such a connection 
is not sure to exist since only the few eigenvalues with largest modulus have a physical 
signification, while we are looking for properties of the entire spectrum. However, our 
numerical results show a connection which we will discuss. We will also give an extension 
of the so-called "disorder variety" to the asymmetric eight-vertex model where the partition 
function can be summed up without Yang-Baxter integrability. We then present an infinite 
discrete symmetry group of the model and an infinite set of algebraic varieties stable under 
this group. Finally, we test all these varieties from the point of view of RMT analysis. 
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This paper is organized as follows: in Sec. we recall the machinery of RMT, and we 
give some details about the numerical methods we use. Sec. |l| is devoted to the eight- 
vertex model. We list the cases where the partition function can be summed up, and give 
some new analytical results concerning the disorder variety and the automorphy group of 
the asymmetric eight-vertex model. The numerical results of the analysis of the spectrum 
of transfer matrices are presented in Sec. [IV]. The last section concludes with a discussion. 



II. NUMERICAL METHODS OF RMT 
A. Unfolding of the Spectrum 

In RMT analysis one considers the spectrum of the (quantum) Hamiltonian, or of the 
transfer matrix, as a collection of numbers, and one looks for some possibly universal sta- 
tistical properties of this collection of numbers. Obviously, the raw spectrum will not have 
any universal properties. For example, Fig. [I] shows schematically three densities of eigen- 
values: for a 2d Hubbard model, for an eight- vertex model and for the Gaussian Orthogonal 
Ensemble. They have clearly nothing in common. To find universal properties, one has to 
perform a kind of renormalization of the spectrum, this is the so-called unfolding operation. 
This amounts to making the local density of eigenvalues equal to unity everywhere in the 
spectrum. In other words, one has to subtract the regular part from the integrated density of 
states and consider only the fluctuations. This can be achieved by different means, however, 
there is no rigorous prescription and the best criterion is the insensitivity of the final result 
to the method employed or to the parameters (for "reasonable" variation). 

Throughout this paper, we call Ei the raw eigenvalues and e$ the corresponding unfolded 
eigenvalues. Thus the requirement is that the local density of the q's is one. We need to 
compute an averaged integrated density of states p(E) from the actual integrated density of 
states: 

p{E) = ^J*J£8{e-E i )de, (2) 

and then we take q = Np(Ei). To compute p(E) from p(E), we have performed a running 
average: we choose some odd integer 2r + 1 of the order of 9-25 and then replace each 
eigenvalue by a local average: 

1 i+r 



E' 



2r + 



J E E 3 , (3) 



and p(E) is approximated by the linear interpolation between the points of coordinates 
(Ehi). We compared the results with other methods: one can replace each delta peak in 
p(E) by a Gaussian with a properly chosen mean square deviation. Another method is to 
discard the low frequency components in a Fourier transform of p(E). A detailed explanation 
and tests of these methods of unfolding are given in Ref . P . Note also that for very peculiar 
spectra it is necessary to break it into parts and to unfold each part separately. Also the 
extremal eigenvalues are discarded since they induce finite size effects. It comes out that of 
the three methods, the running average unfolding is the best suited in the context of transfer 
matrices, and it is also the fastest. 
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B. Symmetries 



For quantum Hamiltonians, it is well known that it is necessary to sort the eigenvalues 
with respect to their quantum numbers, and to compare only eigenvalues of states belong- 
ing to the same quantum numbers. This is due to the fact that eigenstates with different 
symmetries are essentially uncorrelated. The same holds for transfer matrices. In general, a 
transfer matrix T of a classical statistical mechanics lattice model (vertex model) depends on 
several parameters (Boltzmann weights Wj). Due to the lattice symmetries, or to other sym- 
metries (permutation of colors and so on), there exist some operators S acting on the same 
space as the transfer matrix and which are independent of the parameters, commuting with 
T: [T({wi}), S] = 0. It is then possible to find subspaces of T which are also independent of 
the parameters. Projection on these invariant subspaces amounts to block-diagonalizing T 
and to split the unique spectrum of T into the many spectra of each block. The construction 
of the projectors is done with the help of the character table of irreducible representations 
of the symmetry group. Details can be found in [p|,|10|1. 

As we will discuss in the next sections, we always restricted ourselves to symmetric trans- 
fer matrices. Consequently the blocks are also symmetric and there are only real eigenvalues. 
The diagonalization is performed using standard methods of linear algebra (contained in the 
LAPACK library). The construction of the transfer matrix and the determination of its 
symmetries depend on the model and are detailed in Sec. [II [B| for the eight-vertex model. 



C. Quantities Characterizing the Spectrum 

Once the spectrum has been obtained and unfolded, various statistical properties of the 
spectrum are investigated. The simplest one is the distribution P(s) of the spacings s = 
e?+i — Q between two consecutive unfolded eigenvalues. This distribution will be compared 
to an exponential and to the Wigner law ([I]). Usually, a simple visual inspection is sufficient 
to recognize the presence of level repulsion, the main property for non-integrable models. 
However, to quantify the "degree" of level repulsion, it is convenient to use a parameterized 
distribution which interpolates between the Poisson law, the Wigner law. From the many 
possible distributions we have chosen the Brody distribution p], ch. 16.8]: 

Pp(s) = Cl exp (-c 2 / +1 ) (4a) 

with 

and ci = (1 + (3)c 2 . (4b) 

For (3 = 0, this is a simple exponential for the Poisson ensemble, and for (3=1, one recovers 
the Wigner surmise for the GOE. This distribution turns out to be convenient since its 
indefinite integral can be expressed with elementary functions. It has been widely used in 
the literature, except when special distributions were expected as at the metal insulator 
transition Minimizing the quantity: 



c 2 



4 



Jo 



(5) 



yields a value of f3 characterizing the degree of level repulsion of the distribution P(s). We 
have always found <j)(P) small. When —0.1 < (3 < 0.1, the distribution is close to a Poisson 
law, while for 0.5 < f3 < 1.2 the distribution is close to the Wigner surmise. 

If a distribution is found to be close to the Wigner surmise (or the Poisson law), this 
does not mean that the GOE (or the Diagonal Matrices Ensemble) describes correctly the 
spectrum. Therefore it is of interest to compute functions involving higher order correlations 
as for example the spectral rigidity 0: 



A 3 (E) = 




where (. . .) Q denotes an average over the whole spectrum. This quantity measures the 
deviation from equal spacing. For a totally rigid spectrum, as that of the harmonic oscillator, 
one has Af c (E) = 1/12, for an integrable (Poissonian) system one has Ag 01 ^) = E/15, 
while for the Gaussian Orthogonal Ensemble one has A^ OE (E) = -^(log(£ l ) — 0.0687) + 
0(E~ l ). It has been found that the spectral rigidity of quantum spin systems follows 
A^ 01 (E) in the integrable case and A^ OE (E) in the non- integrable case. However, in both 
cases, even though P(s) is in good agreement with RMT, deviations from RMT occur for 
As(E) at some system dependent point E*. This stems from the fact that the rigidity As(E) 
probes correlations beyond nearest neighbours in contrast to P(s). 



III. THE ASYMMETRIC EIGHT- VERTEX MODEL ON A SQUARE LATTICE 

A. Generalities 



We will focus in this section on the asymmetric eight-vertex model on a square lattice. 
We use the standard notations of Ref. 
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The eight-vertex condition specifies that only 
vertices are allowed which have an even number of arrows pointing to the center of the 
vertex. Fig. |2] shows the eight vertices with their corresponding Boltzmann weight. The 
partition function per site depends on these eight homogeneous variables (or equivalently 
seven independent values): 



Z(a, a', b, b' , c, c 



d,d') 



(7) 



It is customary to arrange the eight (homogeneous) Boltzmann weights in a 4 x 4 i?-matrix: 



K 



( a d\ 

b c 

d V 

\d f a' J 



(8) 



The entry IZij is the Boltzmann weight of the vertex defined by the four digits of the binary 
representation of the two indices i and j. The row index corresponds to the east and south 
edges and the column index corresponds to the west and north edges: 
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P 



■ V 



a 



When the Boltzmann weights are unchanged by negating all the four edge values the model is 
said symmetric otherwise it is asymmetric. This should not be confused with the symmetry 
of the transfer matrix. Let us now discuss a general symmetry property of the model. 
A combinatorial argument JR] shows that for any lattice without dangling ends, the two 
parameters c and d can be taken equal, and that, for most regular lattices (including the 
periodic square lattice considered in this work), d and d! can also be taken equal (gauge 
transformation |13[). Specifically, one has: 



b, b' , c, d, d, d!) 



b, b', Vcd, Vcd, V dd', V dd' 



(9) 



We will therefore always take c = d and d — d! in the numerical calculations. In the 
following, when d and d' are not mentioned it is implicitly meant that d = c and d! = d. 
Let us finally recall that the asymmetric eight-vertex model is equivalent to an Ising spin 
model on a square lattice including next nearest neighbor interactions on the diagonals and 
four-spin interactions around a plaquette (IRF model) |OJPl. However, this equivalence is 



not exact on a finite lattice since the L x M plaquettes do not form a basis (to have a cycle 
basis, one must take any LxM-1 plaquettes plus one horizontal and one vertical cycle). 



B. The Row- To-Row Transfer Matrix 

Our aim is to study the full spectrum of the transfer matrix. More specifically, we 
investigate the properties of the row-to-row transfer matrix which corresponds to build up 
a periodic L x M rectangular lattice by adding rows of length L. The transfer matrix Tl is 
a 2 L x 2 L matrix and the partition function becomes: 

Z(a, a', 6, c, d) = Tr [T L (a, a', 6, &', c, d)} M . (10) 

However, there are many other possibilities to build up the lattice, each corresponding to 
another form of transfer matrix: it just has to lead to the same partition function. Other 



widely used examples are diagonal (-to-diagonal) and corner transfer matrices ||12|| . 

The index of the row-to-row transfer matrix enumerates the 2 L possible configurations 
of one row of L vertical bonds. We choose a binary coding: 

L-l 

a = a ^ = K' • • • > a L-l) (ii) 
i=0 

with cti G {0, 1}, corresponding to arrows pointing up or to the right and 1 for the other 
directions. One entry T a ^ thus describes the contribution to the partition function of two 
neighboring rows having the configurations a and (3: 



6 



L-1 



(12) 



With our binary notation, the eight-vertex condition means that w(/ii, a^fy, fii+i) = if 
the sum /Xj + a« + $ + is odd. Therefore, the sum (|12|) reduces to exactly two terms: 
once /io is chosen (two possibilities), /ii is uniquely defined since «o and [3q are fixed and so 
on. For periodic boundary conditions, the entry T a ^ is zero if the sum of all $ and is 
odd. This property naturally splits the transfer matrix into two blocks: entries between row 
configurations with an even number of up arrows and entries between configurations with 
an odd number of up arrows. 

1. Symmetries of the Transfer Matrix 

Let us now discuss various symmetry properties of the transfer matrix. 

(i) When one exchanges the rows a and j3, the vertices of type a, a', b, and V will remain 
unchanged while the vertices of type c and d will exchange into one another. Thus for 
c = d the transfer matrix Ti(a,a' ,b,b' ,c,d) is symmetric. In general the symmetry of the 
row-to-row transfer matrix is satisfied for c = d' and d = d . In terms of the equivalent IRF 
Ising model, condition c = d means that the two diagonal interactions J and J' (confer to 
Ref. [0) are the same: the Ising model is isotropic and therefore its row-to-row transfer 
matrix is symmetric, too. This coincidence is remarkable since the equivalence between the 
asymmetric eight- vertex model and the Ising model is not exact on a finite lattice as already 
mentioned. 

(ii) We now consider the effect of permutations of lattice sites preserving the neighboring 
relations. Denote by 5* a translation operator defined by: 

S\a , «i, . . . , OiL-i) = |«i, • • • , ol-i, ao) • (13) 

Then we have: 

{aS~ l \Ti{a : a, b, b\ c, d)\S(3) = (a|7z(a, a ,b,b' , c, d)\(3) , (14) 

and therefore: 

[T L {a,a',b,b',c,d),S\ = . (15) 
For the reflection operator R defined by: 

R\a 0} ai,..., (Xl-i) = \a L -i, ■ ■ ■ , «i, «o) , (16) 

we have: 

{aR- x \T L {a, a', 6, 6', c,d)\Rfi) = (at\T L (a,a',b,b',d,c)\P) . (17) 
Thus R commutes with T only for the symmetric case c = d: 

[T L (a,a',b,b',c,c),R] = . (18) 
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Combination of the translations S and the reflection R leads to the dihedral group T>l- These 
are all the general lattice symmetries in the square lattice case. The one dimensional nature 
of the group V L reflects the dimensionality of the rows added to the lattice by a multiplication 
by T. This is general : the symmetries of the transfer matrices of d-dimensional lattice 
models are the symmetries of (d — l)-dimensional space. The translational invariance in 
the last space direction has already been exploited with the use of the transfer matrix itself 
leading to Eq. ([II]). 

(iii) Lastly, we look at symmetries due to operations on the dynamic variables themselves. 
There is a priori no continuous symmetry in this model in contrast with the Heisenberg 
quantum chain which has a continuous SU(2) spin symmetry. But one can define an operator 
C returning all arrows: 

C\a , ai,..., = |1 - a , 1 - a 1} . . . , 1 - a L -i) • (19) 

This leads to an exchange of primed and unprimed Boltzmann weights: 

{aC^T^a' ,b,b' ,c,d)\Cj3) = (a\T L (a',a,b',b,c,d)\(3) , (20) 

Thus for the symmetric eight-vertex model (Baxter model) the symmetry operator C com- 
mutes with the transfer matrix: 

[T L {a,a,b,b,c,d),C] = . (21) 



2. Projectors 

Once the symmetries have been identified, it is simple to construct the projectors of 
one row of each irreducible representation of the group T>l (details can be found in |9|,^0[]). 



When L is even, there are four representations of dimension 1 and L/2 — 1 representations 
of dimension 2 (i.e. in all there are L/2 + 3 projectors). When L is odd, there are two one- 
dimensional representations and (L — 1)/2 representations of dimension 2, in all (L — 1)/2 + 2 
projectors. 

For the symmetric model with a = a' and b = b', there is an extra Z% symmetry which 
doubles the number of projectors. 

Using the projectors block diagonalizes the transfer matrix leaving a collection of small 
matrices to diagonalize instead of the large one. For example, for L = 16, the total row-to- 
row transfer matrix has the linear size 2 L = 65536, the projected blocks have linear sizes 
between 906 and 2065 (see also Tabs. | and U). 

As already mentioned, the block projection not only saves computing time for the di- 
agonalization but is necessary to sort the eigenvalues with respect to the symmetry of the 
corresponding eigenstates. In summary, when c = d, the row-to-row transfer matrix is sym- 
metric leading to a real spectrum. Its symmetries have been identified. This is a fortunate 
situation since restriction c = d does neither prevent, nor enforce, Yang-Baxter integrability 
as will be explained in the following section. 
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C. Integrability of the Eight- Vertex Model 

We now summarize the cases where the partition function of the eight-vertex model 
can be analyzed and possibly computed. These are the symmetric eight-vertex model, the 
asymmetric six- vertex model, the free-fermion variety and some "disorder solutions" . 



1. The Symmetric Eight- Vertex Model 

Firstly, in the absence of an 'electrical field', i.e. when a = a', b = b', c = c', and d = d', 
the transfer matrix can be diagonalized using the Bethe ansatz or the Yang-Baxter equations 
This case is called the symmetric eight- vertex model, also called Baxter model 



One finds that two row-to- row transfer matrices T^(a, 6, c, d) and T^(a, 6, c, d) commute if: 

A(a, b, c, d) = A(a, 6, c, d) (22a) 

T{a,b,c,d) = r(d,b,c,d) (22b) 

with: 

ab — cd 

rM ' c ' rf) = ^T^' (23a) 

n 2 _|_ h 2 — r 2 _ A 2 

Note that these necessary conditions are valid for any lattice size L. One also gets the same 
conditions for the column-to-column transfer matrices of this model. Thus the commutation 
relations lead to a foliation of the parameter space in elliptic curves given by the intersection 
of two quadrics Eq. (p3|), that is to an elliptic parameterization (in the so-called principal 
regime fll2"|): 



a = p sn (77 — v) (24a) 

b = p sn (77 + v) (24b) 

c = p sn (277) (24c) 

d = —p k sn (277) sn (77 — v) sn (77 + v) (24d) 

where sn denotes the Jacobian elliptic function and k their modulus. 

It is also well known that the transfer matrix T(a, b, c, d) commutes with the Hamiltonian 



of the anisotropic Heisenberg chain [15 



n = - E + W°3+i + J^Ui ( 25 ) 

i 

if: 

1 : T(a, b, c, d) : A(a, b, c, <f) = J x : J y : J z . (26) 

This means that, given the three coupling constants J x , J y , and J z of a Heisenberg Hamil- 
tonian, there exist infinitly many quadruplets (a, b, c, d) of parameters such that: 
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[T(a, b, c, d),H(J x , J y , J z )\ = 



(27) 



Indeed the three constants J x , J y , and J z determine uniquely r] and k in the elliptic param- 
eterization fl2"I|) and the spectral parameter z/ can take any value, thus defining a continuous 
one-parameter family. Not only T and H, commute for arbitrary values of the parameter v, 
but 7i is also related to the logarithmic derivative of T at v = r\. In this work, we examine 
only regions with the extra condition c = d to ensure that T is symmetric, and thus that 
the spectrum is symmetric. Using the symmetries of the eight-vertex model, one finds that 
the model (a,b,c,d), with c = d mapped into its principal regime, gives a model (a,b,c,d) 
with a — b. In terms of the elliptic parameterization this means sn (rj — v) — sn (rj + v) or 
v = 0. 

In summary, in the continuous one-parameter family of commuting transfer matrices 
T(y) corresponding to a given value of A and T, there are two special values of the spectral 
parameter v: v = 77 is related to the Heisenberg Hamiltonian Ti,(l, T, A), and for v = the 
transfer matrix T{y) is symmetric (up to a gauge transformation). 



2. Six- Vertex Model 



The six-vertex model is a special case of the eight-vertex model: one disallows the two 
last vertex configurations of Fig. 0, this means d = d' = 0. Both, the symmetric and 
asymmetric six-vertex models, have been analyzed using the Bethe ansatz or also the Yang- 



Baxter equations [|12|,|n||I7]] . We did not examine this situation any further since condition 
c = d to have a real spectrum (see paragraph [III B l| (i)) leads to a trivial case. 



3. Free-Fermion Condition 

Another case where the asymmetric eight-vertex model can be solved is the case where 
the Boltzmann weights verify the so-called free-fermion condition: 

aa ' + W = cc' + d<£ (28) 

For condition (p£f) the model reduces to a quantum problem of free fermions and the partition 
function can thus be computed P~8"| , |TP| . 



The free-fermion asymmetric eight- vertex model is Yang-Baxter integrable, however the 
parameterization of the Yang-Baxter equations is more involved compared to the situation 
described in section [HI C 1|: the row-to-row and column-to-column commutations correspond 



to two different foliations of the parameter space in algebraic surfaces. 

It is also known that the asymmetric eight- vertex free-fermion model can be mapped onto 
a checkerboard Ising model. In Appendix we give the correspondence between the vertex 
model and the spin model. The partition function per site of the model can be expressed 
in term of elliptic functions E which are not (due to the complexity of the parameterization 
of the Yang-Baxter equations) straightforwardly related to the two sets of surfaces param- 
eterizing the Yang-Baxter equations or even to the canonical elliptic parameterization of 
the generic (non free-fermion) asymmetric eight- vertex model (see Eqs BBI in the following, 
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see also p5[). The elliptic modulus of these elliptic functions E is given in Appendix [A] 
as a function of the checkerboard Ising variables as well as in the homogeneous Boltzmann 
weights (a, a', b, b', c, c', d, and d!) for the free-fermion asymmetric eight-vertex model. 

Finally, we remark that the restriction c = d is compatible with the condition (|28[) and, in 
contrast with the asymmetric six-vertex model, the asymmetric free-fermion model provides 
a case where the row-to-row transfer matrix of the model is symmetric. 



4- Disorder Solutions 

If the parameters a, a', b, b', c, and d are chosen such that the i?-matrix (|§|) has an 
eigenvector which is a pure tensorial product: 





then the vector: 





(29) 



(30) 



(2L factors) is an eigenvector of the diagonal(-to-diagonal) transfer matrix Tl, usually simply 
called the diagonal transfer matrix. The corresponding eigenvalue is A = X 2L , with 

aa' — bb' + cc' — dd! 
A= (a + a') - (b + V) ■ (31) 

However, the eigenvalue A may, or may not, be the eigenvalue of largest modulus. This 
corresponds to the existence of so-called disorder solutions |20| for which some dimensional 



reduction of the model occurs [21]. Condition ( ]2"P| ) is simple to express, it reads: 
A 2 + B 2 + C 2 + D 2 + 2AB - 2AD - 2BC - 2CD = 

(A + B-C-D)(a + b)(a' + b') -(A- D)(b 2 + b' 2 ) - (B - C)(a 2 + a' 2 ) (32) 

where A = aa', B = bb', C = cc', and D = dd'. Note that in the symmetric case a = a', 
b = b', c = c', and d = d! , Eq. fl3~2|) factorizes as: 



(a — b + d — c) (a — b + d + c)(a — b — d — c) (a — b — d + c) = (33) 

which is the product of terms giving two disorder varieties and two critical varieties of the 
Baxter model. It is known that the symmetric model has four disorder varieties (one of them, 
a + b + c+ d — 0, is not in the physical domain of the parameter space) and four critical 
varieties [12]. The missing varieties can be obtained by replacing TZ by TZ 2 in Eq. (^9|). In 



our numerical calculations we have always found for the asymmetric eight- vertex model that 
A is either the eigenvalue of largest modulus or the eigenvalue of lowest modulus. Finally, 
note that condition (|3^) does not correspond to a solution of the Yang-Baxter equations. 
This can be understood since disorder conditions like ( |3~2"P are not invariant under the action 
of the infinite discrete symmetry group V presented in the next subsection, whereas the 
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solutions of the Yang-Baxter equations are actually invariant under the action of this group 



22,23 



On the other hand, similarly to the Yang-Baxter equations, the "disorder solutions" can 
be seen to correspond to families of commuting diagonal transfer matrices Ti on a subspace 
V of the 2 2L dimensional space on which Tl acts: 

[f L (a, a', b, b', c, d),f L (d, a', b, b', c, d)} \ v = , (34) 

where subscript V means that the commutation is only valid on the subspace V. Actually, 
this subspace is the one-dimensional subspace corresponding to vector fl3~0|). The notion of 
transfer matrices commuting only on a subspace V can clearly have precious consequences 
on the calculation of the eigenvectors and eigenvalues, and hopefully of the partition function 
per site. One sees that the Yang-Baxter integrability and the disorder solution "calculabil- 
ity" are two limiting cases where V respectively corresponds to the entire space where Tl 
acts and to a single vector, namely Eq. (|30|) . 



D. Some Exact Results on the Asymmetric Eight- Vertex Model 



When the Boltzmann weights of the model do not verify any of the conditions of the 
preceding section, the partition function of the model has not yet been calculated. However, 
some analytical results can be stated. Algebraic varieties of the parameter space can be 
presented, which have very special symmetry properties. The intersection of these algebraic 
varieties with critical manifolds of the model are candidates for multicritical points [|24|] . We 
have tested the properties of the spectrum of the transfer matrices on these loci of the phase 
space. 

There actually exists an infinite discrete group of symmetries of the parameter space 
of the asymmetric eight- vertex model (and beyond of the sixteen- vertex model [25]). The 
critical manifolds of the model have to be compatible with this group of symmetries and this 
is also true for any exact property of the model: for instance if the model is Yang-Baxter 

However, 



it is 



integrable, the YBE are compatible with this infinite symmetry group [22 
crucial to recall that this symmetry group is not restricted to the Yang-Baxter integrabil- 
ity. It is a symmetry group of the model beyond the integrable framework and provides for 
instance a canonical elliptic foliation of the parameter space of the model (see the concept 
of quasi-integrability pp]|). The group is generated by simple transformations of the homo- 
geneous parameters of the model: the matrix inversion and some representative geometrical 
symmetries, as for example the geometrical symmetry of the square lattice which amounts 
to a simple exchange of c and d: 



fa d \ 

b c 

c b' 

V d a' J 



(a c\ 

b d 

d V 

\c a' J 



Combining / and t\ yields an infinite discrete group V of symmetries of the parameter space 
This group is isomorphic to the infinite dihedral group (up to a semi-direct product 
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with Z2). An infinite order generator of the non-trivial part of this group is for instance ti -I. 
In the parameter space of the model this generator yields an infinite set of points located 
on elliptic curves. The analysis of the orbits of the group T for the asymmetric eight-vertex 
model yields (a finite set of) elliptic curves given by: 

(aa' + bb' - cd - dd') 2 aa'bb' . . 

const, — — — = const (35) 



aa'bb' ' cc'dd' 



and 



abed 

— = const, — = const, — = const, — = const. 

a' b' d d' 

In the limit of the symmetric eight- vertex model one recovers the well-known elliptic curves 
( P3|) of the Baxter model given by the intersection of two quadrics. Recalling parameteriza- 
tion fl2"4|) one sees that t\ ■ I, the infinite order generator of T, is actually represented as a 
shift by 7] of the spectral parameter: v — * v + rj. 

The group T is generically infinite, however, if some conditions on the parameters hold, 
it degenerates into a finite group. These conditions define algebraic varieties for which 
the model has a high degree of symmetry. The location of multicritical points seems to 
correspond to enhanced symmetries namely to the algebraic varieties where the symmetry 



group T degenerates into a finite group |24j| . Such conditions of commensuration of the shift 
r) with one of the two periods of the elliptic functions occurred many times in the literature 
of theoretical physics (Tutte-Behara numbers, rational values of the central charge and 
of critical exponents |26[]). Furthermore, one can have, from the conformal field theory 
literature, a prejudice of free-fermion parastatistics on these algebraic varieties of enhanced 
symmetry [[27]]. It is thus natural to concentrate on them. We therefore have determined an 
infinite number of these algebraic varieties, which are remarkably codimension-one varieties 
of the parameter space. Their explicit expressions become quickly very large in terms of the 
homogeneous parameters of the asymmetric eight-vertex model, however, their expressions 
are remarkably simple in terms of some algebraic invariants generalizing those of the Baxter 
model, namely: 

J x = Vaa'bb' + Vcc'dd' (36a) 

J y = Vaa'bb' - Vcc'dd' (36b) 

aa' + bb' — cd — dd' 
J z = . (36c) 



Note that, in the symmetric subcase, one recovers Eqs. (P6|). In terms of these well-suited 
homogeneous variables, it is possible to extend the "shift doubling" (rj — > 2rf) and "shift 
tripling" (77 — > 3r]) transformations of the Baxter model to the asymmetric eight-vertex 
model. One gets for the shift doubling transformation: 

(37a) 
(37b) 
(37c) 





- J 2 J 2 


J x J y 


J z ^x 




~ J 2 J 2 


- J 2 J 2 

x °y 


- J 2 J 2 

z °y 


J' z 


- J 2 J 2 


- J 2 J 2 

°z °x 


- J 2 .] 2 

z °y 
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and for the shift tripling transformation: 



J x — ( 2J z J y J x 3J y J z + 2J y J z J x ~\~ ^ x 2 Jy J% Jx Jz Jx J ' Jx (38a) 
Jy = (2J z J y J x — 3J Z J X + J y J x — 2J y J z J x + J y J z + 2 J y J z J X J ■ Jj, (38b) 
^ = + 2JyJ z J x - 3Jpl - 2JpUl + 2J*jy x + Jt4) ■ J z (38c) 

The simplest co dimension-one finite order varieties are: J x = 0, J y = 0, or J z = 0. One 



remarks that J z = is nothing but the free-fermion condition (|28|) which is thus a condition 
for T to be finite. Another simple example is: 

JyJz JxJy JxJz 0? (^^) 

and the relations obtained by all permutations of x, y, and z. Using the two polynomial 
transformations ( |3?| ) and (p8|) one can easily get an infinite number of codimension-one 
algebraic varieties of finite order. The demonstration that the codimension-one algebraic 
varieties built in such a way are actually finite order conditions of T will be given else- 
where. Some low order varieties are given in Appendix 0. In the next section, the lower 
order varieties are tested from the view point of statistical properties of the transfer matrix 
spectrum. 



IV. RESULTS OF THE RMT ANALYSIS 
A. General Remarks 

The phase space of the asymmetric eight-vertex model with the constraint c = d (en- 
suring symmetric transfer matrices and thus real spectra) is a four-dimensional space (five 
homogeneous parameters a, a' , b, b', and c). Many particular algebraic varieties of this four- 
dimensional space have been presented in the previous section and will now be analyzed from 
the random matrix theory point of view. We will present the full distribution of eigenvalue 
spacings and the spectral rigidity at some representative points. Then we will analyze the 
behavior of the eigenvalue spacing distribution along different paths in the four- dimensional 
parameter space. These paths will be defined keeping some Boltzmann weights constant 
and parameterizing the others by a single parameter t. 

We have generated transfer matrices for various linear sizes, up to L = 16 vertices, 
leading to transfer matrices of size up to 65536 x 65536. Tables [I] and [IT] give the dimensions 
of the different invariant subspaces for L = 14 and L = 16. Note that the size of the blocks 
to diagonalize increases exponentially with the linear size L. The behavior in the various 
subspaces is not significantly different. Nevertheless the statistics is better for larger blocks 
since the influence of the boundary of the spectrum and finite size effects are smaller. To 
get better statistics we also have averaged the results of several blocks for the same linear 
size L. 
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B. Near the Symmetric Eight- Vertex Model 



Fig. ^ presents the probability distribution of the eigenvalue spacings for three different 
sets of Boltzmann weights which are listed in Tab. |HI| . Fig. |3|a) corresponds to a point of a 
symmetric eight- vertex model while the other cases (b) and (c) are results for the asymmetric 
eight- vertex model. The data points result from about 4400 eigenvalue spacings coming from 
the ten even subspaces for L = 14 which are listed in Tab. |. For the symmetric model (a), 
using the symmetry under reversal of all arrows, these blocks can once more be splitted 
into two sub-blocks of equal size. The broken lines show the exponential and the Wigner 
distribution as the exact results for the diagonal random matrix ensemble (i.e. independent 
eigenvalues) and the 2x2 GOE matrices. In Fig. |3|a) the data points fit very well an 
exponential, whereas in Figs. |3|b) and |3]c) they are close to the Wigner surmise. In the latter 
cases we have also added the best fitting Brody distribution with the parameter (3 listed in 
Tab. pTT| . The agreement with the Wigner distribution is better for the case (c) where the 
asymmetry expressed by the ratio a /a' is bigger. 

We also have calculated the spectral rigidity to test how accurate is the description of 
spectra of transfer matrices in terms of spectra of mathematical random matrix ensembles. 
We present in Fig. f| the spectral rigidity A^(E) for the same points in parameter space 
corresponding to integrability and to non-integrability as in Fig. The two limiting cases 
corresponding to the Poissonian distributed eigenvalues (solid line) and to GOE distributed 
eigenvalues (dashed line) are also shown. For the integrable point the agreement between 
the numerical data and the expected rigidity is very good. For the non-integrable case the 
departure of the rigidity from the expected behavior appears at E « 2 in case (b) and at 
E ~ 6 in case (c) (in units of eigenvalue spacings), indicating that the RMT analysis is 
only valid at short scales. Such behavior has already been seen in quantum spin systems 
IK?] . We stress that the numerical results concerning the rigidity depend much more on the 
unfolding than the results concerning the spacing distribution. 

Summarizing the results for the eigenvalue spacing distribution and the rigidity, we 
have found very good agreement with the Poissonian ensemble for the symmetric eight- 
vertex model (a), good agreement with the GOE for the asymmetric model (c) and some 
intermediate behavior for the asymmetric eight- vertex model (b). The difference between 
the behavior for the cases (b) and (c) can be explained by the larger asymmetry in case (c): 
case (b) is closer to the integrable symmetric eight-vertex model. 

To study the proximity to the integrable model, we have determined the 'degree' of level 
repulsion (3 by fitting the Brody distribution to the statistics along a path (a = t, a' = 4/i, 

b = b' = 4/5, c = J 5/8) joining the cases (a) and (c) for different lattice sizes. The result is 
shown in Fig. |5|, the details about the number of blocks and eigenvalue spacings used in the 
distributions are listed in Tab. IV. A finite size effect is seen: we always find (3 ~ for the 
symmetric model at a /a' = 1 and increasing the system size leads to a better coincidence 
with the Wigner distribution {(3 = 1) for the non-integrable asymmetric model a ^ a'. So in 
the limit L — > oo we claim to find a delta-peak at the symmetric point a = a' = 2. We also 
have found that the size effects are really controlled by the length L and not by the size of 
the block. However, our finite size analysis is only qualitative. There is an uncertainty on (3 
of about ±0.1. There are two possible sources for this uncertainty. The first one is a simple 
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statistical effect and could be reduced increasing the number of spacings. The second one is 
a more inherent problem due to the empirical parameters in the unfolding procedure. This 
source of errors can not be suppressed increasing the size L. For a quantitative analysis 
of the finite size effects it would be necessary to have a high accuracy on j3 and to vary L 
over a large scale, which is not possible because of the exponential growth of the block sizes 
with L. 

To test a possible extension of the critical variety a = b + c + d outside the symmetric 
region a = a', b = b', we have performed similar calculations along the path (a = t, a' = 4/t, 
b = b' = 4/5, c = 3/5) crossing the symmetric integrable variety at a critical point t = 2. The 
results are the same: we did not find any kind of Poissonian behavior when t ^ 2. We have 
tested one single path and not the whole neighborhood around the Baxter critical variety. 
This would be necessary if one really wants to test a possible relation between Poissonian 
behavior and criticity instead of integrability (both properties being often related). The 
possible relation between Poissonian behavior and criticity will be discussed for spin models 
in the second paper. 

We conclude, from all these calculations, that the analysis of the properties of the un- 
folded spectrum of transfer matrices provides an efficient way to detect integrable models, 
as already known for the energy spectrum of quantum systems [[| ^ . 

C. Case of Poissonian Behavior for the Asymmetric Eight- Vertex Model 

We now investigate the phase space far from the Baxter model. We define paths in the 
asymmetric region which cross the varieties introduced above but which do not cross the 
Baxter model. These paths and their intersection with the different varieties are summarized 
in Tab. [V|. Fig. |6| corresponds to the path (a) (a = 4/5, a' = 5/4, b = b' = t, c = 1.3). 
This defines a path which crosses the free-fermion variety at the solution of Eq. (p8|): t = 
t s = V2.38 and the disorder variety at the two solutions of Eq. (|29|): t = t^ ax « 1.044 
and at t — t™ n ~ 1.0056 (the subscript "di" stands for disorder). See Tab. [V] for the 
intersections with the other varieties. We have numerically found that, at the point t = £J£ ax , 
the eigenvalue (|31~D is the one of largest modulus, whereas at t = t™ n it is the eigenvalue 
of smallest modulus (this is the superscript min or max). The results shown are obtained 
using the representation R = for L — 16 (see Tab. |I1|). After unfolding and discarding 
boundary states we are left with a distribution of about 1100 eigenvalue spacings. One 
clearly sees that, most of the time, j3 is of the order of one, signaling that the spacing 
distribution is very close to the Wigner distribution, except for t = tg « 1.54 and for t close 
to the disorder solutions, where (3 is close to zero. This is not surprising for t = tg since 
the model is Yang-Baxter integrable at this point. The value Pits) is slightly negative: 
this is related to a 'level attraction' already noted in ||. The downward peak is very sharp 
and a good approximation of a 5 peak that we expect for an infinite size. At t = t™ n and 
t = the model is not Yang-Baxter integrable. We cannot numerically resolve between 
these two points. Therefore, we now study paths where these two disorder solutions are 
clearly distinct. For Fig. |8] they are both below the free-fermion point, while for Fig. |7| the 
free-fermion point is between the two disorder solution points. 

In each of the Figs. [7] and [8] are shown the results for two paths which differ only by 
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an exchange of the two weights a «-> a'. In Fig. [7] one clearly sees a peak to j3 slightly 
negative at the free-fermion point at t — 0.8 and another one at one disorder solution 
point t = t^ ax ~ 1.46 for both curves but no peak at the second disorder solution point 
t = £™ n ~ 0.55. It is remarkable that only point t = t™ x yields the eigenvalue of largest 
modulus for the diagonal(-to-diagonal) transfer matrix. Consequently, one has the partition 
function per site of the model at this point. At point t = t™ n , where the partition function 
is not known, we find level repulsion. However, only for path (c) the degree of level repulsion 
(3 is close to unity while for path (b) it saturates at a much smaller value. Another difference 
between the cases (b) and (c) is a minimum in the curve of (3(t) for path (c) at t ~ 1.8 which 
is not seen for path (b). We do not have a theoretical explanation for these phenomena yet: 
these points are not located on any of the varieties presented in this paper. We stress that 
an explanation cannot straightforwardly be found in the framework of the symmetry group 
T presented here since a and a' appear only with the product aa! . It also cannot be found 
in the Yang-Baxter framework, since a and a' are on the same footing in the Yang-Baxter 
equations. 

In Fig. |8| the curves of f3 for the two paths (d) and (e) again coincide very well at the 
free-fermion point at t — ig ~ 1.61. But the behavior is very different for t < 1 where the 
solutions of the disorder variety are located. For the path (d) neither of the two disorder 
points is seen on the curve f3(t) which is almost stationary near a value around 0.6. This 
means that some eigenvalue repulsion occurs, but the entire distribution is not very close to 
the Wigner surmise. On the contrary, for path (e) the spacing distribution is very close to a 
Poissonian distribution ((3(t) ~ 0) when t is between the two disorder points. This suggests 
that the status of eigenvalue spectrum on the disorder variety of the asymmetric eight- vertex 
model is not simple: a more systematic study would help to clarify the situation. 

We now summarize the results from the Figs. [J]-|8|: generally, the statistical properties of 
the transfer matrix spectra of the asymmetric eight-vertex model are close to those of the 
GOE except for some algebraic varieties. We have always a very sharp peak with (3 — > 
at the free-fermion point, often (3 ~ —0.2. All other points with j3 — > are found to be a 
solution of the disorder variety ( |32"1) of the asymmetric eight-vertex model. 

D. Special Algebraic Varieties 

To conclude this section we discuss the special algebraic varieties of the symmetry group 
T. As explained in subsection |U1D| it is possible to construct an infinite number of algebraic 
varieties where the generator is of finite order n: (ti ■ I) n = Id and thus T is finite order. As 
an example, the solutions for n = 6 and n = 16 are given in Appendix [FJ. We have actually 
calculated a third variety, the expression of which is too long to be given (n = 8). We give in 
Tab. [V| the values of the parameter t for which each path crosses each variety tf , tf , t}^ (the 
subscript "fo" stands for finite order and the superscript is the order /). It is easy to verify 
on the different curves that no tendency to Poissonian behavior occurs at these points. We 
therefore give a negative answer to the question of a special status of generic points of the 
algebraic finite order varieties with respect to the statistical properties of the transfer matrix 
spectra. However, one can still imagine that subvarieties of these finite order varieties could 
have Poissonian behavior and be candidates for free parafermions or multicritical points. 
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V. CONCLUSION AND DISCUSSION 



We have found that the entire spectrum of the symmetric row-to-row transfer matrices 
of the eight-vertex model of lattice statistical mechanics is sensitive to the Yang-Baxter 
integrability of the model. The GOE provides a satisfactory description of the spectrum of 
non Yang-Baxter integrable models: the eigenvalue spacing distribution and the spectral 
rigidity up to an energy scale of several spacings are in agreement with the Wigner surmise 
and the rigidity of the GOE matrix spectra. This accounts for "eigenvalue repulsion". In 
contrast, for Yang-Baxter integrable models, the unfolded spectrum has many features of a 
set of independent numbers: the spacing distribution is Poissonian and the rigidity is linear 
over a large energy scale. This accounts for "eigenvalue independence". However, we have 
also given a non Yang-Baxter integrable disorder solution of the asymmetric eight-vertex 
model. For some parts of it the spectrum is clearly Poissonian, too. This suggests that 
the Wignerian nature of the spectrum is not completely controlled by the Yang-Baxter 
integrability alone, but possibly by a more general notion of "calculability" , possibly based 
on the existence of a family of transfer matrices commuting on the same subspace. We have 
also found some "eigenvalue attraction" for some Yang-Baxter integrable models, namely 
for most points of the free-fermion variety. These results could surprise since we do not 
a priori expect properties involving all the 2 L eigenvalues when only the few eigenvalues 
of larger modulus have a physical significance. However, the eigenvalues of small modulus 
control the finite size effects, and it is well known that, for example, the critical behavior 
(critical exponents) can be deduced from the finite size effects. The nature of the eigenvalue 
spacing distribution being an effective criterion, we have also used it to test unsuccessfully 
various special manifolds including the vicinity of the critical variety of the Baxter model. 
We will present in a forthcoming publication a similar study of spin model (rather than 
vertex models). In particular it is interesting to analyze the spectrum on a critical, but 
not Yang-Baxter integrable, algebraic variety of codimension one as it can be found in the 



q = 3 Potts model on a triangular lattice with three-spin interactions [^8]. However, this 
leads models the transfer matrix of which cannot be made symmetric. This will require 
a particular study of complex spectra which is much more complicated. In particular the 
eigenvalue repulsion becomes two-dimensional, and to investigate the eigenvalue spacing 
distribution, one has to analyze the distances between eigenvalues in two dimensions. 
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APPENDIX A: 

We give hereafter the correspondence between the asymmetric eight-vertex model on the 
free-fermion variety and the checkerboard Ising model. The vertex model is specified by 
eight homogeneous parameters. The gauge invariance @, together with the free-fermion 
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b 



condition (p8[) , leads to only four independent parameters. The checkerboard model is 
specified by the four independent usual low-temperature variables X{ = e i . 

a = - xx 2 x 3 x 4 H , a = 1 (Al) 

2 V X1X2X3X4/ 2 \X2X4 X1X3J 

1 / XiXj + x 2 x 3 \ fe/ _ 1 / XiX 2 + x 4 x 3 \ 

2 \x 2 x 3 X1X4J ' 2 VX4X3 :cix 2 / 

1 (X X X 2 X Z X 4 \ ,1 x 2 \ /AoN 

2 V x 4 xix 2 x 3 / 2 V x 2 X1X4X3/ 

^ _ 1 /X4X3X2 Xi \ j _ 1 / X1X4X2 X 3 \ ^ 

2 V a;i £4X3X2/ ' 2 \ x 3 Xi£4£2/ 

The modulus &; c heck of the checkerboard Ising model or equivalently of the asymmetric eight- 
vertex model reads : 

^check = ~i (A5) 
K D 

with : 

= 16 (a^a^a^ + x 4 ) x circular permutations 
k D = (xlxlxlxl — x\x\ — x\x\ + x\x\ + X3X4 + x\x\ — x\x\ — 1) x (circular permutations) 

or equivalently: 

2 16 c'd'cd 

free = Wc'd'cd - 8a'bab' - 2a' 2 b 2 - 2a 2 b' 2 - 2a 2 b 2 - 2b' 2 a' 2 + {a 2 - a' 2 ) 2 + [b 2 - b' 2 ) 2 ( ' 



APPENDIX B: 

We give in this appendix some algebraic varieties where the group T becomes finite. Two 
varieties of order six ((ij - J) 6 = Id): 

cd 3 ±c 2 d 2 + dc 3 - cdbb' 

a± = 7a X liJj ^ B1 ) 

aide ± a bb' 

A variety of order sixteen ((t\ - J) 16 = Id): 

(u 2 - v 2 ){u - v) + uv{u 2 + v 2 ) =0 (B2) 



with: 



cc'dd' (aa' + W - c 2 - d 2 ) 2 _ 1 _ fj x J z - J y J z x 1 
aa'bb — cc'dd') 2 \ J^J,, 



" = )Zui ' -1=1 — -1 (B3) 



aa'bb'iaa' + bb' - c 2 - d 2 ) 2 ( J X J Z + J V J. x 2 



(aa'bb — cc'dd') 2 \ J x J y 



- 1 = 2 - 1 (B4) 



19 



REFERENCES 



* E-mail: hmeyer@crtbt.polycnrs-gre.fr, dauriac@crtbt.polycnrs-gre.fr 
t E-mail: maillard@lpthe.jussieu.fr 
[1] E. P. Wigner, Ann. Math. 53, 36 (1953). 

[2] M. L. Mehta, Random Matrices, 2nd ed. (Academic Press, San Diego, 1991). 
[3] D. Poilblanc, T. Ziman, J. Bellisard, F. Mila, and G. Montambaux, Europhys. Lett. 22, 
537 (1993). 

[4] T. C. Hsu and J.-C. Angles d'Auriac, Phys. Rev. B 47, 14291 (1993). 

[5] G. Montambaux, D. Poilblanc, J. Bellisard, and C. Sire, Phys. Rev. Lett. 70, 497 (1993). 

[6] P. van Ede van der Pals and P. Gaspard, Phys. Rev. E 49, 79 (1994). 

[7] H. Bruus and J.-C. Angles d'Auriac, Europhys. Lett. 35, 321 (1996). 

[8] H. Meyer, J.-C. Angles d'Auriac, and H. Bruus, J. Phys. A: Math. Gen. L (1996). 

[9] H. Bruus and J.-C. Angles d'Auriac, CRTBT preprint (1996). 
[10] H. Meyer, Ph.D. thesis, Univ. J. Fourier, Grenoble, France, 1996. 
[11] I. Varga, E. Hofstetter, M. Schreiber, and J. Pipek, Phys. Rev. B 52, 7783 (1995). 
[12] R. Baxter, Exactly Solved Models in Statistical Mechanics (Academic Press, New York, 
1982). 

[13] A. Gaaf and J. Hijmans, Physica A80, 149 (1975). 

[14] P. W. Kasteleyn, in Fundamental Problems in Statistical Mechanics III, Proc. of the 
1994 Wageningen Summer School, edited by E. Cohen (North- Holland, Amsterdam, 
1975), pp. 103-155. 

[15] B. Sutherland, J. Math. Phys. 11, 3183 (1970). 

[16] E. H. Lieb and F. Y. Wu, in Phase Transitions and Critical Phenomena, edited by C. 

Domb and M. Green (Academic Press, New York, 1972), Vol. 1, pp. 331-490. 
[17] I. M. Nolden, J. Stat. Phys. 67, 155 (1992). 
[18] C. Fan and F. Y. Wu, Phys. Rev. B 2, 723 (1970). 
[19] B. U. Felderhof, Physica 66, 509 (1973). 
[20] M. T. Jaekel and J.-M. Maillard, J. Phys. A 18, 1229 (1985). 

[21] A. Georges, D. Hansel, P. Le Doussal, and J.-M. Maillard, J. Phys. A20, 5299 (1987). 
[22] M. P. Bellon, J.-M. Maillard, and C.-M. Viallet, Phys. Lett. B 260, 87 (1991). 
[23] J.-M. Maillard, J. Math. Phys. 27, 2776 (1986). 

[24] H. Meyer, J.-C. Angles d'Auriac, J.-M. Maillard, and G. Rollet, Physica A 209, 223 
(1994). 

[25] M. Bellon, J.-M. Maillard, and C.-M. Viallet, Phys. Lett. B 281, 315 (1992). 
[26] J.-M. Maillard and R. Rammal, J. Phys. A16, 353 (1983). 

[27] See for example S. Dasmahapatra, R. Dedem, T. R. Klassen, B. M. McCoy, and E. 

Melzer, in Proc. of the Conf. Yang-Baxter Equations in Paris, edited by J.-M. Maillard 

(World Scientific, Singapore, 1993). 
[28] F. Y. Wu and R. K. P. Zia, J. Phy. A 14, 721 (1981). 



20 



TABLES 



R 


k 


A 


Ir 


. even 
a R 










1 


1 


362 


325 


1 


7T 


-1 


1 


288 


325 


2 





-1 


1 


234 


261 


3 


7T 


1 


1 


288 


261 


4 


2vr/7 


* 


2 


594 


585 


5 


4vr/7 


* 


2 


594 


585 


6 


6vr/7 


* 


2 


594 


585 


7 


tt/7 


* 


2 


576 


585 


8 


3vr/7 


* 


2 


576 


585 


9 


5vr/7 


* 


2 


576 


585 



TABLE I. The dimensions clr and degeneracies Ir of the invariant subspaces for L = 14. R is 
an arbitrary label of the representations of the dihedral group, exp(ifc) and A are the eigenvalues 
of the corresponding translation and reflection operators (* means that the corresponding repre- 
sentation is not stable under the action of the reflection operator). The two numbers a^ en and 
a R dd correspond to the blocks between configurations with an even or odd number of up arrows 
(see text). 



R 


k 


A 


Ir 


„even 
a R 


„odd 
a R 








1 


1 


1162 


1088 


1 


7T 


-1 


1 


1033 


1088 


2 





-1 


1 


906 


960 


3 


7T 


1 


1 


1033 


960 


4 


tt/2 


* 


2 


2065 


2048 


5 


vr/4 


* 


2 


2062 


2048 


6 


3vr/4 


* 


2 


2062 


2048 


7 


tt/8 


* 


2 


2032 


2048 


8 


7vr/8 


* 


2 


2032 


2048 


9 


5vr/8 


* 


2 


2032 


2048 


10 


3vr/8 


* 


2 


2032 


2048 



TABLE II. The parameters of the invariant subspaces as in the preceding table for L = 16. 



case a a' b = b' c = d (3 



a) 


2 


2 


0.8 


75/8? 


d 0.790569 


0.02 


b) 


1.8 


2.2222 


0.8 


^5/8 s 


d 0.790569 


0.64 


c) 


1.6 


2.5 


0.8 


\/578 s 


a 0.790569 


0.85 



TABLE III. The Boltzmann weights for the three cases of figures |3j and ^. The last column 
indicates the parameter (3 of the Brody distribution shown in Fig. ||. 
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L dimension aj of the blocks ^(spacings) 

12 122, 176, 174, 166, 165,165 870 

13 190, 315, 315, 315 1000 

14 362, 596, 596, 576 2000 

15 612, 1092, 1094 2600 

16 1162, 1092, 2065 4000 



TABLE IV. Numbers of representations and their dimensions used for the distributions leading 
to Fig. |. 



(a) (b) (c) (d) (e) 



a 
a' 
b 
b' 

c = d 


0.8 
1.25 

t 
t 

1.3 


0.4 
2 

1.5 
1 


2 

0.4 

t 

1.5 
1 


1.2 
0.8 
t 

1.5 
1.3 


0.8 
1.2 

t 

1.5 
1.3 


ts 

j. min 
r di 
j. max 
l di 

E fo 

E fo 

r fo 

r fo 


1.5427 
1.0057 
1.0443 

1.3 
1.5991 
1.6000 


0.8 
0.5520 
1.4608 
0.6667 
0.8141 
0.8141 
0.5237 


0.8 
0.5520 
1.4608 
0.6667 
0.8141 
0.8141 
0.5237 


1.6133 
0.6156 
0.2424 
1.1267 
1.7474 
1.7001 
0.6337 


1.6133 
0.6156 
0.2424 
1.1267 
1.7474 
1.7001 
0.6337 


TABLE V. 


The parameterization of the paths 


(a)-(e) and their intersection points 


with the 



different algebraic varieties (see text). 
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FIG. 1. Schematic examples for the den- 
sity of eigenvalues for the 2d Hubbard, GOE 
random matrices, and transfer matrices of the 
eight-vertex model. 



a a' b b' c c' d d! 
FIG. 2. The Boltzmann weights of the eight 
vertices. 
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FIG. 3. The distribution P{s) of eigenvalue 
spacings for the three sets of Boltzmann weights 



enumerated in Tab. ID. The continuous lines 



give the exponential for a Poissonian spectrum 
and the Wigner distribution (Q) for the GOE, the 
broken lines in figures (b) and (c) are the fitted 



Brody distribution with the j3 given in Tab. ID 
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FIG. 4. The spectral rigidity A^(E) for the 
three sets of Boltzmann weights enumerated in 



Tab. |HI| The lines give the limiting case for a 
Poissonian spectrum and for the GOE. 
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FIG. 5. The degree of level repulsion (3 ob- 
tained by a fit of the eigenvalue spacing dis- 
tribution P(s) to the Brody distribution for 
lattice sizes L = 12, . . . , 16. The Boltzmann 
weights vary along the path (a = t, a' = 4/t, 
6 = 6' = 4/5, c = d = y/5/8) which gives a sym- 
metric model for t = 2. The number of spacings 



used for each size is given in Tab. IV 
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FIG. 6. The degree of level repulsion (5 for 
path (a) (see Tab. |V]) (L = 16 and 1000 spacings) 
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FIG. 7. The degree of level repulsion for 
paths (b) and (c) (see Tab. 0) (L = 16 and 1000 
spacings for (b) and 3000 for (c)) 
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FIG. 8. The degree of level repulsion (3 for 
paths (d) and (e) (see Tab. 0) (L = 14 and 
about 2000 spacings) 



25 



